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SECTION  1 
INTRODUCTION 


The  interaction  of  blast  waves  from  two  or  more  simultaneous 
explosions  can  lead  to  much  higher  pressures  than  would  be  found  from 
the  individual  blast  waves.  This  phenomenon  has  been  called  cumulation 
by  Kuhl  (Reference  1)  in  a  thorough  analysis  of  the  interaction  for  a 
six-burst  situation  described  below. 

A  field  test  designated  MISERS  BLUFF  II-2  (MB  11-2)  was  con¬ 
ducted  by  the  Defense  Nuclear  Agency  at  the  Planet  Ranch  test  site  near 
Lake  Havasu,  Arizona,  30  August  1978,  for  measuring  the  interaction  of 
six  blast  waves  from  the  simultaneous  detonation  of  six  120-ton  ANFO 
charges  arranged  in  a  l-.exagonal  pattern  with  100-meter  separation. 

Reference  1  ptv?sents  a  comprehensive  analysis  of  the  basic 
interactions  involved  in  this  six-burst  blast  interaction.  Reference  1 
also  develops  a  semi-empirical  method  for  prediction  of  the  peak  shock 
overpressures  ana  makes  comparisons  with  the  measured  data  from  the 
MB  11-2  test.  The  need  for  a  hydrocode  to  be  used  in  prediction  of  the 
interaction,  particularly  to  obtain  overpressure-time  waveforms,  is 
noted  in  Reference  1. 

The  objective  of  the  present  investigation  has  been  to  deter¬ 
mine  whether  the  REFLECT  hydrocode  would  have  the  potential,  if  modified, 
for  computing  the  interaction  of  blast  waves  from  multiple  nuclear 
bursts  simultaneously  detonated  at  ground  level  and  to  indicate  the 
modifications  needed. 

A  principal  feature  of  the  REFLECT  code  is  the  representation 
of  the  refracted  (reflected)  shock  waves  in  a  two-burst  interaction  as 
"sharp"  shocks  using  the  Rankine-Hugoniot  (R-H)  equations.  If  blast 
waves  from  additional  bursts  interact  with  these  blast  waves,  as  for  the 
six-burst  array  being  studied,  it  is  expected  that  the  additional  shock 
waves  might  also  have  to  be  represented  in  such  a  calculation  by  the 
R-H  equations  in  order  to  obtain  satisfactory  accuracy. 

The  approach  taken  has  been  to  make  a  limited  modification  to 
the  REFLECT  code  so  that  a  preliminary  comparison  could  be  made  of  the 
code  results  with  the  MB  II-2  test  data.  Only  the  two  shock  waves  that 
are  currently  modelled  in  the  code  as  R-H  shocks,  the  first  refracted 
shock  and  the  first  Mach  shock,  would  be  modelled.  The  shock  waves 
appearing  later  would  then  be  smeared  shocks  in  these  calculations. 

This  comparison  would  enable  an  assessment  to  be  made  of  the  potential 
of  the  REFLECT  code  were  one  or  more  of  the  additional  shock  waves  to  be 
modelled  as  R-H  shocks. 
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Sharp  R-H  shock  representation  is  expected  to  be  very  important 
to  accurate  calculation  of  blast  wave  interaction  from  multiple  bursts. 
Without  R-H  shock  modelling,  the  shocks  become  smeared  as  compressions 
over  a  number  of  cells  introducing  strong  numerical  viscous  effects  into 
the  computed  interaction.  Several  other  improvements  to  the  REFLECT 
code  should  also  be  considered  for  the  calculations,  and  they  will  be 
discussed  and  summarized  below. 


SECTION  2 

REFLECT-3  HYDROCODE 


2.-1  REFLECT-1  MODEL  FOR  TWO-BURST  INTERACTION. 

The  blast  wave  pattern  for  six  charges  in  a  hexagon^’!  pattern 
is  shown  in  Figure  1  for  a  time  shortly  after  the  waves  begj  Inter¬ 

act.  The  blast  waves  all  first  interact  at  the  charge  bisec. 

Because  the  cnarges  are  equally  spaced  and  equidistant  from  the  array 
center,  the  wave  pattern  is  symmetrical  on  the  two  sides  of  ach  of  the 
six  charge  bisectors. 

The  REFLECT-1  code  (Referenci-  2)  computes  the  flow  field 
within  the  shaded  region  for  a  plane  of  symmetry  which  is  in  this  case 
the  charge  bisector  plane.  The  solution  is  axisymmetric  about  the  line 
joining  the  two  charges.  A  sketch  illustrating  the  REFLECT-1  solution 
and  its  relationship  to  the  six-burst  problem  is  shown  in  Figure  2.  The 
code  computes  the  flow  field  between  the  refracted  blast  shock  and  the 
charge  bisector  plane  until  the  time  that  the  refracted  shock  Intercepts 
the  plane  of  the  charge  radial. 

The  REFLECT  code  employs  the  Godunov  differencing  scheme 
(Reference  3)  with  an  expanding  cell  system.  The  cell  system,  Figure  2, 
is  axisymmetric  about  the  axis  passing  through  two  adjacent  burst 
centers.  Each  cell  is  a  torus  centered  about  the  axis.  The  cells  are 
arranged  in  columns  with  the  generators  for  each  column  emanating  from  a 
common  point  on  the  axis  on  the  opposite  side  of  the  charge  bisector 
plane. 

The  refracted  blast  shock  Rj^  is  one  boundary  of  the  cell 
system.  The  charge  bisector  plane  is  the  second  boundary  of  the  cell 
system.  After  the  Mach  shock  forms  between  the  intersection  of  the 
shocks  Ip  and  Rp  and  the  bisector,  it  becomes  the  third  boundary  of  the 
cell  system.  Further  description  of  the  cel.l  system  is  given  in  Reference  2. 

The  Godunov  method  conserves  mass,  momentum  and  energy  in  each 
cell.  The  equation  of  state  for  air  is  taken  from  Brode  (Reference  4). 

The  REFLECT-1  code  then  computes  the  flow  within  the  region  in 
Figure  2  bounded  by  the  refracted  shock  Rp,  the  charge  bisector  plane 
and  the  Mach  shock  Mp  after  it  forms.  Figure  3.  The  boundary  condition.s 
for  the  calculation  are  the  following.  At  the  sliock  Rp  the  properties 
and  Clow  velocities  of  the  blast  wave  immediately  outside  of  th.e  shock 
provide  tlie  boundary  conditions.  At  the  bisector  plane  the  impermeable 
surface  conditions  apply.  At  the  Mach  shock  Mp,  when  it  forms,  the 
atmospheric  conditions  apply  at  the  boundary. 

The  cell  layout  used  in  the  present  calculations  is  the  same 
layout  employed  in  Reference  2.  There  are  24  columns  of  cells  between 
the  axis  and  the  triangular  comer  cell  at  the  juncture  of  the  Rp  wave 


CHARGE  BISECTOR 


Figure  2.  REFL.ECT-l  model  of  two-burst  blast  interaction. 
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with  the  charge  bisector.  Tlie  number  of  cells  in  a  column  decreases 
stepwise  from  tl>e  axis  to  the  corner  cell.  Also,  the  number  of  cells 
within  a  column  increases  with  time  as  the  cell  system  expands.  The 
maximum  number  of  cells  in  a  column  is  limited  to  lb. 

After  the  Mach  sliock  is  formed  the  triangular  corner  cell 
is  replaced  by  one  or  more  four-sided  cells.  This  cell  array  is  caJled 
the  outer  mesh  system  in  Reference  2.  The  maximum  number  of  cells 
reached  in  the  outer  mesh  system  in  the  prc "ent  calculations  was  16  in  a 
4x4  array. 

The  number  of  cells  used  in  the  present  calculations  was 
limited  to  the  capability  of  the  REFLECT  code  developed  in  Reference  2. 

A  recommended  improvement  to  be  emphasized  below  would  be  to  increase 
the  number  of  cells  significantly. 

2-2  SIX-BURST  INTERACTION. 

Up  to  the  i^oint  in  the  calculation  where  the  shock  Ri  readies 
the  charge  radial  tlie  REFLECT-3  code  is  the  same  as  the  REFLECT-1  code 
of  Reference  2,  with  the  exception  of  improvements  made  to  the  code 
regarding  iteration  methods,  tolerances  and  other  associated  techniques, 

When  the  Rj^  shock  reaches  the  charge  radial  the  changes  made 
for  the  REFLECT-3  code  apply.  The  cell  layout  for  this  period  is  shown 
in  Figure  3,  where  each  cell  shown  in  the  sketch  represents  four  colls 
in  the  calculation. 

The  adaptation  made  to  the  REFLECT-1  code  for  this  problem  has 
been  to  limit  the  expansion  of  the  cell  system  to  the  boundary  provided 
by  the  charge  radial.  Therefore  each  column  of  cell  cross-sections 
shown  in  Figure  2  expands  until  it  reaches  the  cdiargc  radial,  as  sliown 
in  Figure  3.  In  the  inner  nine  columns  shown  in  the  sketch  the  cells 
are  bounded  by  the  charge  radial.  The  outer  four  columns  are  bounded  by 
the  reflected  shock  R][^  and  the  Mach  shock  As  the  incident  shock  I] 

moves  closer  to  the  array  center,  the  remaining  columns  one  by  one  reach 
the  charge  radial.  During  the  whole  period,  of  course,  the  width  of  the 
columns  increases . 

When  the  cells  reach  the  charge  radial,  a  pressure  wave 
reflects  back  towards  the  charge  bisector.  'J'his  is  the  refracted  shock 
wave  from  the  next  clockwise  burst,  90  degrees  from  the  bisector. 

The  sequence,  of  waves  from  the  time  the  blast  waves  first 
interact  is  sho\^m  in  Figure  4.  The  waves  are  sketched  in  this  figure 
essentially  in  the  shape  obtained  in  the  REFLECT-3  results.  The  shock 
waves  that  are  computed  correctly  as  R-H  shocks  are  identified  by  double 
lines  and  the  other  shocks  by  single  lines. 
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Adjacent  blast  waves  first  interact  on  “he  charge  bisectors  at 
284.1  feet  from  the  array  center.  Figure  4a.  The  shock  pattern  after 
the  initial  interaction  is  shown  in  Figure  4b.  Up  to  this  point  all  the 
shock  waves  are  computed  by  the  REFLECT-3  code  as  R-H  shocks. 

The  pattern  after  the  Mach  shock  forms  is  shown  in  Figure  4c. 
The  refracted  shock  R2  is  not  represented  as  an  R-H  shock  ir  the  REFLECT~3 
code,  so  it  smears  over  several  cells.  The  results  of  the  REFLECT-3 
calculations  indicate  the  triple  point  intercepts  the  charge 

radial  at  24.5  feet  from  the  array  center.  This  compares  well  with  the 
24  feet  given  by  the  method  of  Kuhl  (Reference  1) . 

After  the  triple  point  l;i^-Rj^-Mj  intercepts  the  charge  radial, 
the  Mach  shock  M2  is  formed  with  the  triple  point  Mi-R2“M2.  Because  the 
shock  R2  is  smeared  in  the  present  form  of  the  REFLECT-3  code,  the 
triple  point  becomes  spread  out  along  the  incident  shock  ML-M2.  The 
smeared  triple  point  therefore  would  be  expected  to  roach  the  charge 
bisector  ahead  of  the  time  it  would  if  instead  it  were  formed  by  an 
R^  shock  that  is  sharp. 

The  arrival  of  the  principal  part  of  the  smeared  triple  point 
was  identified  from  tlie  REFLECT-1  results  by  the  abrupt  rise  in  pressure 
behind  the  M^  sliock  as  being  at  9  feet.  This  compares  well  with  the 
arrival  of  the  triple  point  given  by  Reference  1  an  8  feet. 

The  shock  R2  refracts  across  the  charge  bisector  as  the  R3 
wave  shown  in  Figure  4e  and  the  Much  shock  M-j  is  formed.  The  shock  R3 
is  smeared  in  the  REFLECT-3  calculations,  but  tlie  front  of  the  implod’jng 
wave  cotitinuos  to  be  formed  by  the  R-ll  shock  in  the  code.  This  imploding 
sliock  has  the  triple  point  M2-R3-M2  according  to  the  model  of  Reference  1. 

The  triple  point  M2-R3-M3  is  smeared  in  tlie  calculations 
because  the  shock  R3  is  smeared.  The  arrival  of  the  leading  edge  of  the 
smeared  mone  at  the  charge  radial  is  indicated  in  the  computed  results 
by  a  marked  rise  in  pressure  behind  shock  M2  at  the  racial  at  3  feet 
from  the  array  center.  From  there  inward  the  results  iudicatud  that  the 
pressure  on  tlie  rad  il  behind  the  incident  shock  iiicrennod  to  esHeiitinlly 
equal  the  pressure  at  the  bisector. 

The  REFLECT-3  calculations  were  continued  until  the  incident 
shock  on  the  bisector  reached  within  0.43  feet  of  the  array  center.  At 
that  point  the  cell  system  was  frozen.  From  then  on  the  waves  moved 
through  the  mesh  with  all  the  waves  being  smeared.  At  this  point  the 
cell  system  had  four  cells  along  the  imploding  shock  between  the  charge 
bisector  and  charge  radial.  Each  cell  had  a  0.056-ft  front  on  the 
imploding  shock  and  a  length  in  the  radial  direction  of  0.20  feet. 

The  calculations  were  continued  until  0.073  seconds  after  the 
scaled  AFWL  l-KT-STD(REV)  burst.  This  corresponds  uo  about  0.080 
seconds  for  the  120-ton  ANFO  bursts. 
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Thn  PFAB  code  has  been  developed  for  computing  the  shock  wave 
that  radiates  from  the  array  center  after  the  imploding  wave  r  laches  the 
center.  (This  wave  is  called  t;he  main  reflected  wave  in  Reference  1). 
This  code  is  described  in  the  Appendix.  The  code  computes  the  region 
inside  of  the  outward-facing  shock  and  the  sector  bounded  by  the  charge 
bisector  and  radial.  The  boundary  conditions  for  this  zone  arc  the 
conditions  ahead  of  the  outward-facing  shock,  given  by  che  REFLECT-3 
calculations,  and  the  impermeable  surface  conditions  at  the  charge 
bisector  and  radial.  The  PFAB  code  for  computing  this  region  has  been 
prepared  and  debugged. 

A  computer  routine  is  required  to  input  the  necessary  data 
from  the  REFLECT-3  code  into  the  PFAB  code.  The  REFRA  mating  code  was 
developed  for  use  with  the  REFLECT- 1  code  output,  and  minor  modifications 
are  required  for  its  use  with  tlie  REFLEGT-3  output. 

There  is  u  point  concerning  tlie  three-dimensionality  of  the 
REFLECT-3  model  that  sliould  be  noted.  The  vertical  plane  from  the 
ground  along  the  charge  radial  is  represented  in  tlic  REFLECT-3  code  as  a 
slightly  curved  surface.  This  is  due  to  the  fact  that  the  cell  geometry 
of  the  REFLECT-3  code  is  axiuymmetric  about  the  axis  joining  tlie  charges. 
In  the  code  the  boundary  for  the  charge  radial  is  a  conical  surface 
having  its  apex  at  the  charge.  The  charge  radial  on  the  ground  is  a 
generatrix  of  tlie  cone.  The  includiid  angle  of  the  cone  is  120  degrees, 
so  the  radius  of  curvature  uormai  to  the  surface  is  very  largo  compared 
with  the  distance  from  the  radial  to  the  bisector,  and  the  ratio  of  the 
two  becomes  infinite  approaching  the  array  <’euter.  Because  the  radius 
of  curvature  becomes  increasingly  largo  relative  to  the  distance  between 
the  charge  radial  and  uhurge  bisector  as  the  incident  shock  approaches 
the  charge  center,  calculations  have  indicated  tliut  the  curvature  is  not 
expected  to  have  a  significant  influence  on  the  results,  particularly 
approaching  the  array  center. 
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SECTION  3 
RESULTS 


The  test  results  of  the  MISERS  EhUFF  II-2  test  are  presented 
in  Figures  5  to  12.  These  results  have  been  photocopied  from  Reference  1. 

Tlie  free-uJr  blast  data  are  an  Jnpiit  to  the  KKFLECT--3  code. 

'J^o  blast  models  are  compared  witli  the  Lest  data  In  Figures  5  and  6. 

The  records  of  the  overpressures  measured  in  the  MB  11-2  test 
are  presented  in  Figures  6  and  7  for  stations  along  the  charge  radials 
and  bisectors,  respectively.  The  records  have  been  reproduced  from 
Reference  1.  The  stations  are  Identified  by  the  radial  distance  from 
the  array  center  in  meters  and  angic  from  one  of  the  charge  bi.sectors  in 
degrees;  e.g,,  Station  25-30  is  25  meters  from  the  array  center  and 
30  degrees  from  the  baseline  charge  bisector. 

The  output  of  the  REFl.ECT-3  calculations  is  compared  wltli 
the  test  results  in  Figures  6  and  7  fur  tlie  ciiargo.  radials  and  bisectors, 
respectively.  The  times  for  the  KEFLECT--3  resjults  presented  in  tliese 
two  figures  have  been  adjusted  to  match  the  measured  arrival  of  tlwj 
first  sliock  of  each  lu'essurc  record.  The  incident  shock  pressure-range 
rosults  are  compared  with  the  test  data  in  Figure  8.  The  pressure-range 
results  for  tlio  main  reflected  wave  are  compared  with  the  test  data  in 
Figures  9  and  10.  The  sliock  positions  as  a  function  of  time  are  compared 
with  the  test  rosults  in  Figures  11  and  12. 

3-1  ULAbT-WAVii  MODEL. 

3-J. .1  Shock  Overpressure. 

Tlie  REFLECT  code  uses  tlie  data  from  the  undisturbed  blast  wave 
from  u  single  burst  an  input  to  tlie  code  calcul.itions .  The  culcuiatijiis 
are  as  accurate  as  this  input. 

Tlic  incident  sliock  overpressure,  Ap^,  given  by  two  methods  Ls 
compared  in  Figure  5  with  data  from  the  MB  fI-2  test  results.,  The 
MB  lf-2  data  used  here  for  comparisons  have  been  obtained  from  Reference  1. 
The  abscissa  is  the  distance  from  the  center  of  tlie  6-bur, at  array. 

The  crosses  in  Figure  5  represent  the  ANFO  eaiculutiun  liy  AFWI. 
for  the  MB  11-2  test  as  gJven  in  Reference  5.  The  shock  overpressure 
presented  here  for  these  data  was  determined  by  the  cnrii'ut  method 
recommended  by  Reforene??  6  wherein  the  shock  overpressure  is  taken  as 
the  peak  computed  pressure,  and  the  location  of  the  shock  is  taken  as 
the  point  whore  the  computed  pressure  is  u  .e-half  of  the  peak  pressure. 
Previously  Reference  7  had  reconuiiended  a  4-point  extrapolation  of  the 
prassures  to  obtain  the  shock  overpressure.  That  method  generally  gives 
a  shock  overpressure  about  2  to  A  percent  higher  .-.han  for  the  curve 
shown.  At  the  82-ft  station  the  ANFO  calculation  give,s  a  shock  over¬ 
pressure  that  is  10  percent  higher  than  a  mean  of  the  data.  For  the 
41-ft  station  the  ANFO  calculation  shows  good  agreement  with  the  test 
data . 
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The  second  blast  model  was  obtained  from  scaling  the  AFWL 
l-KT-STD(REV)  of  Reference  8  for  a  free-air  yield  of  0.67  KT  for  the 
ambient  conditions  of  the  MB  II-2  test.  The  shock  overpressure  for 
this  model  is  also  compared  in  Figure  5  with  the  MB  11-2  test  results. 

The  agreement  with  the  test  data  is  very  good  at  both  tlie  82  and  41-ft 
stations . 

3-1.2  Overpressure  Wave  Form. 

At  test  stations  on  the  radials  at  a  distance  of  82  feet  from 
the  array  center.  Figure  6a,  the  overpressure  during  approximately  the 
first  0.025  seconds  following  the  arrival  of  the  shock  is  produced  by 
the  blast  wave  from  the  charge  located  on  the  radial.  For  this  period 
the  overpressure  wave  forms  for  the  ANFO  calculations  and  the  scaled 
1-KT-STD (REV)  models  can  be  compared  against  the  test  data.  Two  records 
have  been  used  in  the  comparison:  from  Stations  25-30  (30°  radial)  and 
25-210  (210°  radial).  For  Station  25-30  the  ANFO  computed  values,  using 
at  the  shock  the  4-point  extrapolation,  are  about  10  percent  higher  (as 
a  fraction  of  the  shock  overpressure)  than  the  test  data  during  about 
the  first  10  milliseconds,  they  essentially  agree  with  the  test  data 
during  the  middle  of  the  period  and  fall  about  5  percent  below  toward 
the  end.  The  scaled  1-KT-STD (REV)  values  are  in  good  agreement  for 
the  whole  period. 

For  Station  25-210  the  ANFO  calculations  agree  fairly  well 
with  the  test  data.  The  scaled  1-KT-STD (REV)  values  are  about  10  percent 
lower  than  the  test  data  at  early  times  and  about  10  percent  higher  when 
the  period  ends  (with  the  arrival  of  the  Rj^  wave).  It  is  concluded 
from  these  comparisons  that  the  waveform  measurements  at  82  feet  from 
the  array  center  support  the  ANFO  and  l-KT-STD(REV)  blast  models  about 
equally  well. 

Two  stations  were  located  on  charge  radials  at  41  feet  from 
the  array  center:  Stations  12.5-30  and  12.5-150.  The  records  are 
shown  in  Figure  6b.  The  waveform  lasts  only  about  0.005  and  0.003 
seconds,  respectively,  at  the  two  stations.  The  record  for  Station 
12.5-150  indicates  the  Ip  waveform  is  terminated  by  a  double  peak. 

Except  for  that  the  scaled  l-KT-STD(REV)  values  agree  well  with  these 
waveform  data  to  within  the  resolution  of  the  graphs  (about  0.2  to 
0.4-in.  amplitudes). 

These  comparisons  on  the  radials  provide  an  evaluation  of 
the  modelling  of  single  ANFO  blast  waves  at  the  ground  level.  From 
these  comparisons  it  is  concluded  that  the  scaled  l-KT-STD(REV)  blast 
model  agrees  with  the  test  data  within  the  data  scatter  at  82  and  41  feet 
of  the  array  center  until  0.06  seconds.  There  are  no  single-burst  blast 
data  nearer  to  the  array  center  than  the  41-ft  station  or  later  than 
about  0.06  seconds,  because  of  the  blast  interactions  from  the  other 
charges. 
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3-1.3  Selection  oJ;  Scaled  AFWL  1-KT-STD (REV)  Blast  Model. 

Because  of  the  differences  in  shock  overpressure  at  the 
82‘'£t  station,  it  is  concluded  that  the  scaled  AFWL  l-KT-STD(REV)  blast 
model  is  in  some  'hat  better  agreement  with  the  test  data  than  is  the 
ANFO  calculatlonal  results.  As  far  as  the  overpressure  wave  form 
comparisons,  however,  the  two  models  agree  about  equally  well  with  the 
test  results. 


There  is  another  problem  with  using  the  ANFO  calculatlonal 
results.  The  properties  in  the  region  of  the  smeared  shock  must  be 
defined  and  that  would  require  some  exploratory  work. 

Tliereforc  the  decision  was  made  to  use  the  scaled  AFWL 
l-K.T-f!Tl)(REV)  blast  model  witli  a  free-air  yield  of  0.67  KT  as  the  ii.put 
to  the  RKl''LliCT-3  calculations  to  represent  the  120-ton  ANFO  charges  used 
in  the  MB  11-2  test.  It  is  concluded  that  the  ANFO  calculatlonal  results 
support  this  hla.st  model  for  use  in  this  inner  100-ft  region  of  the 
blast  array  witlilu  about  10  percent  in  shock  overpressure. 

3-2  WAVES. 

The  initial  interaction  of  the  blast  waves  occurs  along  the 
hiiiectoru .  The  Rj  waves  are  produced  by  this  interaction.  A  record  for 
thiu  interaction  is  shown  in  Figure  7a  for  Stations  86.6-60  and  86.6-240. 

The  initial  shock  overpressure  for  an  86.6-m  bisector  station 
can  hu  I'aiculated  directly  from  the  proi>ert:ies  of  the  incident  shocks  by 
uslu)',  tlic  R-11  equations.  This  is  e.Msentially  an  exact  calculation, 
dependent  ui'on  the  Htrength  of  the  incident  shocks.  Because  the  two 
blast  models  (ANFO  computed  and  scaled  1-KT-.STU(KEV) )  were  shown  to 
agree  witliin  about  10  I'ercont  witli  the  measurements  of  the  shock 
overijressure  and  overpressure  waveform,  this  calculation  might  be  con¬ 
sidered  for  use  as  an  evaluation  of  tlie  test  simulation.  The  ANFO 
culculati'd  blast  wave  gives  an  initial  shock  overpressure  that  is 
7  percent  higher  than  the.  scaled  1-KT-STD (REV)  value  at  the  86.6-m  station. 

Tlie  test  data  from  Station  86.6-60  lias  a  shock  overpressure 
that  is  2.S  iiorceut  below  tlie  scaled  l-KT-.STU (REV)  value  and  the  test 
data  of  Station  86.6-240-AB  has  a  shock  oveiprossure  that  is  36  percent 
below  that  vaJuc.  Therefore  it  is  concluded  that  the  shock  overpressures 
measured  at  Initial  blast  intercept  appear  lo  be  low  by  about  25  percent 
or  more. 


'I'hese  differences  cannot  be  used  as  a  definite  evaluation  of 
the  experiment  because  of  questions  regarding  use  of  the  blast  models  at 
Llu'  l()4-fl  (5()-m)  blast  range.  Tlic  comiiarlHons  made  in  Figure  5  are 
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for  2^6-^t  and  287-ft  blast  ranges,  respectively,  for  the  82-ft  and 
41-ft  stations.  The  ANFO  calculation  showed  a  shock  overpressure  that 
tends  to  be  high  at  the  lower  blast  radius  (82-ft  station).  The  scaled 
AFWL  1-KT-STD (REV)  shock  overpressures,  although  they  agreed  at  the  246 
and  287-ft  blast  ranges,  are  much  larger  in  the  proximity  of  the  burst 
center  than  for  a  high  explosive.  Therefore  a  definite  conclusion 
cannot  be  reached  at  this  point  on  the  test  accuracy.  But  the  com¬ 
parison  does  point  up  the  possibility  that  the  Lest  overpressure  might 
be  low  by  as  much  as  25  percent  or  more  compared  with  an  ideal  test. 

Tlie  predicted  overpressure  waveform  for  Station  86.6-60  from 
the  REFLECT-3  calculation  are  shown  in  Figure  7a  to  be  within  10  percent 
of  the  test  data  by  one  millisecond  after  shock  arrival.  The  differences 
with  the  test  data  for  the  86.6-240  Station  are  much  greater  until  about 
five  milliseconds  after  the  intercept. 

T'he  record  for  the  164-ft  station  on  the  charge  bisector. 

Figure  7b,  appears  to  indicate  that  one  blast  shock  arrives  about  1  1/2 
milliseconds  ahead  of  the  other,  tending  to  reduce  the  peak  pressure. 

The  REFLECT-3  shock  value  is  29  percent  higher  than  the  measured  peak. 

At  the  82-ft  stations  op  the  charge  bisector.  Figure  7c, 

Station  25-0  has  a  large  shock  Jump  at  0.0448  seconds  followed  by  a 
second  jump  about  0.003  seconds  later.  Station  25-240  shows  a  similar 
sequence,  although  the  overpressure  in  the  later  case  falls  off  about 
15  percent  before  the  second  jump.  The  REFLECT-3  prediction  essentially 
fairs  through  the  test  data  of  Station  25-0  for  about  0.013  seconds 
after  the  arrival  of  the  wave.  The  second  jump  in  the  overpressure 
may  be  due  to  an  irregular  type  of  Mach  reflection.  Irregular  Mach 
reflectior  resuJ.ts  in  a  compressive  peak  following  the  Mach  shock  due 
tc  the  interaction  of  a  slipline  with  the  (virtual)  reflecting  surface. 

At  about  0.058  seconds  the  REFLECT-!  overpressures  in  Figure  7c 
begin  to  climb  above  the  test  values.  This  departure  is  attributed  to 
the  early  arrival  of  the  leading  edge  of  the  smeared  R2  shock.  The 
R3  wave  arrives  early  as  a  result.  It  is  expected  that  this  effect 
could  be  eliminated  by  making  the  R2  shock  a  sharp  R-H  shock  in  the 
REFLECT-3  code,  as  are  the  Rj^  and  shocks.  In  general,  however,  the 
predictions  shown  in  Figure  7c  are  in  fair  agreement  with  the  test  record. 

For  the  52-ft  station  on  the  charge  bisector,  Fi.gure  7d,  the 
REFLECT-!  prediction  agrees  well  with  the  first  peak  of  the  trace.  The 
measured  pressure  trace  has  a  slight  dip  until  about  0.06  seconds,  but 
the  agreement  with  the  REFLECT-3  prediction  otherwise  is  good.  Beginning 
at  0.06  seconds  the  predictions  again  show  the  effect  of  the  early 
arrival  of  the  smeared  R2  shock  followed  immediately  by  the  R3  wave. 

At  41  feet  on  the  charge  bisector.  Figure  7e,  the  REFLECT-3 
prediction  agrees  well  with  the  first  peak  measured  value  for  the 
wave  for  Station  12.5-240  shown  in  the  lower  trace.  However,  the  measured 
peak  for  Station  12.5-0  is  about  20  percent  lower  and  has  a  flat  waveform 
for  about  3  milliseconds,  followed  by  a  sudden  dip,  instead  of  the  more 
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classical  wave  form  of  Station  12.5-240.  At  about  0.060  seconds  tlie 
overpressure  predicted  by  the  REFLECT-3  code  rises  above  tVie  measured 
values  for  both  the  0-deg  and  240-deg  stations  due  to  the  smeared 
R2-R^  wave. 

The  comparison  for  the  20-ft  stations  on  the  240-deg  bisector 
is  shown  in  Figure  7f.  The  REFLECT-3  predictions  for  the  wave  are 
25  to  50  percent  higlier  than  the  measured  values  for  Station  6-240.  How¬ 
ever  the  predicted  values  fall  below  the  measurements  of  Station  6-180 
by  up  to  25  percent.  So  the  REFLECT-T  prediction  falls  between  tlie  measure¬ 
ments.  Several  factors  make  it  appear  that  the  calibration  of  the  record 
for  Station  6-180  may  be  in  error;  this  will  be  discussed  in  Section  3-5. 

The  predicted  overpressures  begin  to  rise  at  about  0.065  seconds  due  to 
the  early  arrival  of  the  smeared  R2-R3  wave. 

3-3 

The  Rj  wave  in  the  REFLECT-3  code  reaches  the  radial  as  a 
sharp  shock  as  shown  in  Figure  4c.  It  refracts  across  the  radial  as 
the  smeared  R2  wave  in  the  present  form  of  the  code.  The  R]  shock 
appears  to  refract  at  the  radial  as  a  regular  shock  reflection.  The 
R2  wave  following  the  intercept  of  the  R-]  shock  is  indicated  on  the 
traces  in  Figure  6. 


In  Figure  6a  the  R2  wave  comparisons  are  shown  for  82  feet  from 
the  array  center.  The  arrival  time  of  tlie  R1-R2  wave  for  the  REFJdiCT-3 
prediction  is  within  0.0008  seconds  of  the  measured  arrival  time  for 
Station  25-30  (note:  the  measured  variation  in  arrival  time  between  the 
30  and  210-dQg  test  stations  is  0.0004  seconds).  This  is  good  agreement. 
A  second  shock  arrives  0,00085  and  0.0021  seconds  after  the  first  K] 
shock  at  the  two  82-ft  radial  stations.  The  REFLECT-3  results  also  show 
a  rise  following  the  shock  arrival  to  a  second  peak  about  0.004  seconds 
later.  In  fact  both  of  the  test  traces  show  a  third  peak  a  few  milli¬ 
seconds  behind  the  second.  The  second  peak  may  indicate  a  delayed 
arrival  of  one  of  the  R|^  waves  at  eacli  of  the  stations,  and  tlie  third 
peak  may  indicate  an  irregular  type  of  interaction  of  tlie  .sliocks. 

The  REFLECT-3  results  in  Figure  6a  again  fair  through  the  test 
data  for  the  R2  wave  fairly  well  up  to  about  0.056  seconds.  At  that 
point  the  REFLECT-3  values  become  higlier,  which  is  attributed  again  to 
early  arrival  of  later  waves  due  to  shock  smearing.  Resolution  of  the 
sharply  peaked  R2  shock  wave  siiown  in  Figure  6a  would  also  lequire 
reduction  in  the  size  of  the  cells  for  the  REFLECT-3  code. 

Figure  6b  shows  the  comparison  at  the  41-ft  distance.  For 
Sration  12.5-30  the  predicted  arrival  of  the  Ro'Rf  wave  is  0.0010 
seconds  after  the  measured  arrival,  which  is  within  the  0.00125-s 
difference  measured  between  the  arrivals  for  the  two  stations.  Again 
the  REFLECT-3  results  fair  through  the  test  results  until  0.067  seconds. 
The  predictions  then  rise,  which  is  again  attributed  to  the  wave  smearing 
by  the  code,  in  this  case  by  the  smeared  main  reflected  (M-R)  wave. 
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TUe  in-(.'dicLi.oM!i  aj^uiu  do  noL  sliow  Llio  alwir|)  i)oak  ol  llu'  tual.  ilntn,  bill 
that  would  reiiuiro  tho  amulloi'  coIIh.  TIio  prudiotloiiH  do  not  a}'i‘co 
quite  aa  well  with  tho  tost  data  for  Station  12,5-150  but  that  station 
has  an  anomalous  wavol'orm  wJ.Lh  tlio  two  sinallor  waves  proeedint',  the 
apparent  H1-K2  wave  rIvIur  a  sip.nif  Leantly  difforont  overall  waveform 
than  for  Station  12.5-50. 

At  the  20-ft  stations  on  the  ehurgo  radlals,  l'’ij»uro  be,  tho 
llliFLECT-5  overpressure  predietions  have  a  waveform  somewliat  similar  to 
that  of  Station  b  dl.O  for  about  0,004  seconds  after  arrival  of  tlie  shock 
M2.  At  Station  b-2/0  the  M2  shock  is  preceded  by  a  weaker  shock,  whereas 
at  the  150  (not  shown)  and  210-dui.i  radlals  tlie  M2  shock  arrives  first. 

The  level  of  the  predictions  is  about  25  to  50  percent  hiyliur  than  the 
measured  values,  l)ut  there  are  dil'ferences  of  +20  percent  between  tlie 
measured  U2  overpressure  wave  forms  for  tlie  three  stations  at  tliis 
range,  so  it  is  not  clear  tliat  the  measured  interactions  at  this  dis¬ 
tance  from  tlie  array  center  are  accurate  enough  for  much  better  resolution, 
but  it  does  appear  that  tho  trend  of  the  predicted  values  is  higlnr  than 
for  the  test  data  at  this  range. 

At  about  O.Obb  seconds  the  RKFIiKCT-T  iiredietions  in  Figure  be 
begin  to  rise  sii'iiif icantly  above  the  test  values.  Tills  I'ollows  the 
trend  at  tho  more  outward  stations  attributed  to  the  efl'ects  of  shock 
smearing  and  cell  size. 

3-4  R2-U3  WAVE. 

Thu  H2  wave  intercepts  the  charge  bisector  and  reilec.ls 
(regular)  as  an  K3  shock,  as  shown  in  Figure  4e,  both  shocks  are  smeared 
in  the  REFl.ECT-3  code.  The  R3  wave  appears  at  the  b2-it  stations  shown 
in  Figure  7c,  but  it  occurs  too  lute,  for  complete  eomparison  with  the 
test  data.  Comparison  of  the  REi''hE(lT-3  results  with  the  Rq  wave  form  in 
Figures  7d  to  g  show  that  tho  REFLECT-3  wave  is  smeared  and  would  need 
a  sharp  K2  shuck  and  smaller  cells  to  reproduce  tlie  relatively  sharp 
peak  of  the  R3  wave. 

Figure  4e  shows  that  stations  on  charge  biseetors  within  y 
feet  of  tho  array  center  would  be  expected  to  have  the  Mach  shock  M3 
arrive  first.  Carefui  Insiiectioii  of  tlie  trace.  Cor  Station  1.22-lHU 
shows  a  steep  initial  shock  ;lump  very  nearly  equal  to  the  REFhEOT-l 
shock  .Itiiiip  indteated  by  the  first  point.  Tlie  HEFhEC'l'-3  values  continue 
to  climb  us  do  the  test  data,  but  the  REFhECT-3  values  eontinue  to  climb 
due  to  the  early  arrival  of  the  smeared  M-U  wave.  The  initial  shock 
jump  is  not  quite  us  distiiiet  for  the  truce  of  Station  1.22-240,  l)ut  the 
REFLECT-3  results  compare  similarly. 

3-5  M-R  WAVE. 

The  blast  waves  from  charges  that  are  located  more  than  90 
degrees  away  around  the  hexagon  produce  the  wave  identified  here  as  the 
M-R  wave,  taken  again  from  the  terminology  in  Reference  1  where  it  is 
called  the  main  reflected  wave.  For  points  on  the  radlals,  as  shown  in 
Figure  6,  the  M-R  wave  is  produced  by  three  charges.  For  points  on  the. 
bisectors,  as  shown  in  Figure  7,  it  is  produced  by  two  charges. 
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AL  Lliu  Htutii)i\n  within  20  icet  uf  Lho  ufi'uy  cuntor  the  M-H 
wave  ut'i.'ivoH  oatly  oiiounh  to  allow  compnrlHona  ol.'  the  KliM.,li;CT“3  I'ofjultw. 
Thuao  compar iBonti  at'o  uhown  in  Fip.ureM  6c  and  7f  to  i. 

Tlio  M-K  wave  ptodlctod  by  the  K1'!1''IJ';CT“3  code  1h  of  coiirao 
Him'ai'c'd,  but,  wiLli  one  excoption,  the  nuiftnitude  of  the  ptuik  ovin'prenauroM 
Hcnoral. ly  apiteeH  faily  well  wltli  tlie  teut  ni'dulLi).  The  exception  in  for 
Station  6-lHU  which  ia  believed  to  luive  an  erroneonn  c/\lib>.‘ntion ,  which 
la  Muppoi'ted  by  several  other  f(5atureH  (the  rapid  de.c.reaHe  li\  the  over- 
preaimre  of  tlie  Nj  ahoek  alonp,  tlu».  26U“tle|',  cliart'c  l)inec.tor  that  would 
renuit  between  20  and  10  feet  of  tlie  array  center  and  by  c.omparlMon  of 
reouitn  with  tlione  lor  Station  6-260). 

•3“0  ARRAY  OKNTRR. 

Tlu!  eoiiiparlfum  for  near  the  center  of  the  array  iw  whowu  in 
KlPiio'e  '/!.  The  RI'H'TdiOT-’l  reaultn  are  for  a  iiolnt  on  the  blt-iector  0.63 
1  rout  tho  array  center.  Thu  M.^  iihoek  Jump  predicted  by  the  RERLECT-3 
code  1m  V/(I  |im1  ('i.Oi  iiiei'apaMca  1  a) .  it  compnroM  well  wit.h  the  mhock 
jump  iiieaMured  at  the  two  iitatiomi,  an  indicated  by  the  firut  ahurp  riae, 

The  peak  overpruunure  pted feted  by  the  RKVld'XlT-O  codiii  at 
0.6')  ieel  111  1()J')  pul  (li.l  niotpipanca  1  ii)  ,  which  In  7i  percent  of  flio 
peakn  ol  ahout  22(iO  piii  (Ib.ti  nieipipat.ca  1  n)  iiK'atiui'ed  I  roiii  iliu  roeordii  ol 
the  two  Htaliunn  in  Fippire  7i  at  tho  array  center.  Thiu  fa  very  Rood 
aRreumeut  in  view  ul  the,  ri'lattvely  larne  tii/u  of  the  celln  employed  in 
the  RMl'T.KCT-eode  and  l  lu'  nmuariup,  of  the  R2  and  R'^  nhoek  wavea,  iTie 
wave  I'unii  of  the  Rlil'TiliOT-3  ovurpruanure  la  relatively  aimilar  to  the 
meamired  wave  form  of  tlwi  M-K  wave  exeeiiL  for  the  Hrnuter  width  aiul 
reduced  I'eak.  The  eoiiiiiariuon  in  very  uncourayiiia  fur  the  pol.iiuitia'i  of 
the  code . 

3-7  INOiUliNT  SHOOK  OVliRRRIiSSURlSS . 

The  shock  overpreaunres  for  the  incident  uliock  at  |Hrlnt;a  on 
radials  and  bisectors  within  100  feet  ol'  the  array  eentnr  are  plotted  in 
flRuru  S  reproduced  from  Rufurenee  1.  The  Rlil‘T.l!0T-3  valucii)  are  pltjttod 
for  compar  I  noil . 

The  Rl';i''hl'',(;']'-')  radial  values  from  100  feet  Inward  to  26  feet 
are  the  scaled  i-KT-S'fU(RHV)  input  values  to  the  code,  so  they  do  not 
provide  a  coiiiparisuii  for  code  evaluation.  On  the  b  face  tor  at  b2  teet 
the  RKKhliOT-3  values  are  about  10  to  23  ))ore,eiu:  hi  ('her  than  given  for 
the  teat  data.  However  the  teat  data,  Ki.Ruro  7e,  have  a  seeotul  peak 
occurr inR  only  about  0.003  ueeundn  after  the  incident  rihoek,  and  tho 
Klil'’LK0T-3  value  ia  only  7  percent  hlgiier  than  l.lie  Hecoiid-peak  value  for 
Station  2,3-0'*'Al). 

At  32  feet  on  the  biauetor  tlie  KliFl.l''iCT-3  values  ai'/ree  v;itb  the 
test  peak  value.  At  61  feet  on  the  bisector  the  RKFhECJT-l  ebook  over¬ 
pressure  aRrecB  with  the  test  value  for  tlie  260-deR  station.  'J'he  shock 
value,  for  the  O-deg  station  Is  lower  by  16  pereoi'iU,  but  Lhltt  record, 
f'lRurc  7e,  has  a  peculiar  waveform,  so  tho  value  is  questioned. 
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Al;  20  locL  on  llu?  l).Un.'('Luf  lliu  pvoil  1  cit'd  vijJno  ol  1A(^  pH  I 
raiiw  bolwcui'i  Llu!  iiicaMui'od  valuoa  a)'  90-*  10/  pAJ.  and  Idd^-'/IO  pal  lor  Llic 
2A0  and  IHO-dcp,  1)  laoi'./oi'ii ,  it'HiHu;livttly.  liowcvor,  llic  cal  Ibrat  li'U  ol 
Lbc  laLluL  vocord  upiioarn  quoatiounbU’,  aa  dJaciiHHcd  In  dcclJoii  .'l-b,  an 
lUu  p)  od U',1  ion  may  Iju  nlinhily  hlp.h.  Tlio  KKFl.''’.i)T-.l  curve  IVr  Llia 
lilaotLi)!'  in  I'i p.ur'.'  I!  liau  a  aiuali  clan  l)ctwni'u  Llu'  22  diul  JH-l'i  sLallona. 
'lliia  vl.ac  In  ai i), IhnLcd  to  Ma'  advance  arrival  ol  l.lio  amnared  Jl| -,H2“l'l2 
ta'tpiu  poiuL.  Ai  10  lout,  on  Lti>*  biacc  Cic  (.lu>  pvcilicLcil  valuca  arc  ,‘iliown 
in  ,1/1)', arc  (i  Lo  bn  aboiil  ,’U/  I'/crccot  liinbui  tliar  (die  mcai'orcd  pe.ilt,  valuna 
Irom  Llic  liSO  and  240-dOH  biacoiora. 

Al.  a  (Jlataucc  ol  alnitii.  9  loot.  I'rom  ilic  array  iianLcr  ihc 
Rl!d''l.iHOT'-.‘j  incidoivL  aliock  ovorprcaaiiro  cm  Lho  Ijiac.cLoi  lioi'Jaa  to  i;Jliiib< 
Thin  cmrrt'Hp'Uida  I'lilli  i;lu.'  jiroiicucu  ol  (.ho  M'j  wave  abown  In  l'M,^ii);'c  /ic« 
Ik'tauiiio  I'hu  1<2  and  K;)  w/ivok  arc  uoi  iuodc.i.ind  in  l>.)i’.l''Ll.'10T-.’j  al.  Lhl.s  Lime 
.‘in  !(■  II  nliuckn,  t.lio  protliciocl  Inc.ldont  lyhook  ovorpi'ciuniro  i.iuirn.iacn 
amoolhiy  Inal.tiad  nl  witb  (.lie  ahoelt  Jump  tlvAi  would  bo  o’l'eci.od  in  iliu 
acl.ual  caac,  Al  h  I'oni  Irom  l;.liu  ar),7iy  cmilior  i.lic  pri'dicLcd  iiuvidunt 
aliock  ovo vpr>iunu/'c  la  about  iwo'kbirthi  ol  tlw.'  vi/lnen  i‘op(.irted  in 
Ru ioi'itmcw  1.  lluwovor,  the  rt'ctud  lor  Sbation  l-22**i90  iian  a  ,'vtcup 
iiiltla.!  vlae  Lhai  oorveapondn  wo,ii  wlili  tho  Ui'U''l,.i''.0T'“.'i  value,  and  thia 
r.iau  may  bo  tho  ticttniL  M;(  ahock,  and  tin.*  va  iuu  rei>''i'tcil  ;ln  Uciornnci.!  ,1 
iiwy  biii  thill  o.i  a  poiik  iilUu'  the  Mj  Mlmclt  an'.ival, 

Tho  HKKhKOT-.'i  iihi,iiik  ovovi’i’dimufo  on  the  r.ul.laj,  lor  tho  liic.I.doiiL 
ahock  t'Ofu’.hea  .’iy:)  pw,!,  wlniin  l.ho  aliock  in  0./i',)  loot  Irom  tliii  array  eontev. 
Al,  ihut  poinl;  tho  oo,il  ayiitem  waa  Iro/tra.  On  tho  hlauctor  tho  Inoldont 
ahock  waa  wiLliiii  0.41  loot  ol  Lho  array  caiiti.ir  with  an  ovoi'prcaiuire 
tyl  veil  hy  l.ho  KliFhi'!0T“‘,'j  cotUi  ol  H/O  pai. 

It  ia  coiicluiUid  ihai  nho  lacldiiiit  alun'k  atrcmptlia  on  the 
l)i.tuu:Lura  are  i.'rudic,Led  well  into  41  tcuU,.  l'’rom  20  to  10  loot  tho 
predict, iuna  appear  Lo  run  aomowhat  lii)(h  by  aa  imic.li  aa  about  30  ptri'ooi'iL , 
which  may  btt  duo,  to  early  arrJvai.  ol  Llic  R.)  wave  duo  lo  wmeariiip,.  At 
4  loot  tho  Rl','l''hIi'.(;T'“3  I'unulta  .ippo.arn  to  .UKri.'O  wo.ll  with  Lho  data.  Tlie 
nuiiio  ap.ruoi'iuiit  aI.ao  cxLnLn  near  tliw  arr/Jiy  center.  Jio,  the.  RJi.l'hhOT'-.'i 
reHu'itrv  lor  tho  lueidonit  ahock  uverprownuro  on  tiu'  radial  wiipoara  to 
a)';i"wo  wel.l  wl.Lli  what  are  boli 'wed  to  be  tlici  boat  data,  except  whuro 
amuar,tn)/,  ol  tho  wave  n|i)u'.uirfj  (ni  produce  Inlghor  predicted  ovorprenaiirea 

Alouy  the.  radial  the  Ui';l''Id!;0T- ,3  Incident  .ahock  ovorproaaure  la 
lU'Gn  111  Flijuiro  il  to  tncroauo  at  about  the  24-lt  atatlon  whore  Llu* 
wave,  nliuwn  Indlcatoiii  a  Jump.  Ayyiiln,  the  incident,  ahock  In  tho  ).U'',l'’hl':0T-3 
model  if)  u  U-'H  ahock  but  the  relle.cLod  K2  wave  ia  dllluae,  no  the  rise 
would  be  expected  to  be  not  ao  auddcii  aa  in  the  tent.  At  22  1/2  leot 
l.lie  RKF,i,,.'.';CT--3  value.  .ii)truen  with  the  /.thock  interact, ton  analy.Mi«  predict. ion 
ol  Relerence  1.  Al:  the  2,0-lt  fU-ntlon  the  KE1''LKCT~.3  Incident  .vhock 
overprenuuro  ia  about  2.3  to  '.U)  percout.  higher  tlmu  the  iiiuaBured  peak, 
vicluca.  However,  an  moiitionocl  in  Section  .3'-3  there  are  aignif .leant 
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il  1 1  k’ t't'iu'i’H  Jn  Lhu  wnvDlunun  muaHiifocl  on  lIr>  lIuh’o  radlaJM  (150,  210  (iiid 
270  do)’,i‘0(.!H ,  I'l'fijUru  (k').  In  I'/U’ t ,  Lho  wnvul'oi'm  lor  Lho  21U-do)’  rndliil 
oonLlnuoft  Lo  olliiili  alLor  ahock  arrival,  roaohlnr,  n  poak  only  O.OOIB 
tinooiula  alLor  altock  arrival  (lial  aj’rooa  wall  wl.Lli  tlio  prodloLcd  ahock 
ovu  i.'proHHuru ,  On  Lhu  270"do}j  radial  IL  aiipo/ira  Lho  M2  wavu  liaa  noL 
ovortakun  Llio  liu'l.dunL  aliock  at  Llio  20-t'L  aLatlon. 

Tlio  RMFLI'XIT'-'l  prodlcLod  InaldonL  alionk  ovarpri'ftHuri'  on  I  ha 
radial,  Mippira  H,  conLiiuia.u  It)  Inai-atiaa  j>railna  1.1  y  Lo  about  Ilia  l-lL 
al.al.  ion.  At  Lha.  lO-lL  and  4-l‘L  ataLionn  it  in  aiiouL  .SO  to  11. 'j  paraant 
liipjiar  Lliaii  ilia  maai-inrad  [latdt  vuluaa,  buL  a.omparlatma  vkLli  ilia  LaiiL  wava 
I'oriiiii  liava  not.  baaii  iiiada.  Al  A  laaL  tlia  iiiaaauratl  nlioak  ovarpraHHiiraa 
ranp,a  i  1.M111  102  Lo  '322  pal  (Ualarana.a  I)  and  Llia  UMMliKtlT-'J  valtia  in 
'.10)3  pal.  At  till,'  !l-“l  L  Mtatlon  Llia  Ml';i'’i.MflT-'l  vahia  raaaliaa  anaanllally 
Llia  aaiiia  valuo  an  on  Llia  blaaator  and  raiiuiina  naarly  Llia  aaina  Inward  Lo 
U.A'J  laaL  Ironi  Lha  array  aaiiLur  wliara  tlia  iiiaali  waa  aLoppad. 

IL  ia  aoii(;l.udad  Llial.  Lha  inaldaiiL  ahoak  aLraiip.Lli  ia  pratllaLad 
wall  whara  Lha  ahoa.k  llraL  I'oiiiia  aL  abouL  Lha  27t-rL  atalJoii.  Tlia 
Kl';i''l,l''.LlT-’l  valua  conl.inuau  Lo  a  I  l.iiih  wharaaa  Lha  LaaL  valuaa  daaraaaa 
Ini'.ard  Lu  tha  lO-l'i  ataLlon,  but  Lha  tl  1 1  raranaaa  ara  iinraao  1  vad .  I'miii 
A  t'uuL  Inward  Lha  pradla.Lioaa  ara  mipportad  b>  noma  of  Lha  taiiL  raaiilLa, 
biiL  Lha  LaaL  n.'nulLH  dllt\r  a ip,n I  (  1  a.aiiL ly . 

■J-K  MAI  N  RKl'T.KCTKll  BlIDOK  OVKRl’HMUSUHliS . 

.Shock  ovarprouuuraH  pradlatad  for  tha  main  I'al'iaaLad  wava  by 
Lha  Rl''d'T,K(IT“'j  auda  arc  cumiiaratl  wll.h  Lhu  MU  11-2  data  In  I'innran  '■j 
and  10.  For  Lhu  a.huri;a  radial  (hara  ia  omi  poliil.  lor  coinpar  1  aon  which 
iu  aL  20  laaL.  Tha  ap/raania.iL  wiLh  Lha  Lu.cl;  raanlLa  lu  axa.allaiiL. 

For  Lha  tdiatj'o  biuoa.Lor,  Fip.nra  10,  thara  ara  prudlcL ioiia  lor 
Ihraa  loauLioiia  i.o  .10  laaL.  Tha  a p,).' Daman t  aL  10  and  A  luat  la  a.xaallaiiL. 
AL  U.4J  luat;  tha  pratllaLion  la  aboiiL  2'1  paraant.  low.  'll.  la  bwlluvad 
l.lu.iL  a  aniaLlur  i.all.  alisa  wcaila  Impvovu  Lha  prail  la.Lloii ,  iiartlanlarl.y 
Lua.auuu  ul  Lhu  aiaapiiana  of  lha  paak  (FIgnra  71). 

:j-y  llUOOM  I’kONT  TRAJMtlTOUlliS. 

Tha  pradioLad  li'ajaa.ft.ry  ol  Lha  InalduiU.  ahoak  along  a  aliarga 
radial  li»  a.oniparad  wiLli  Lha  taut  daia  In  Figure  11.  Tnw.'ird  l.o  tha  23- 
IL  lUallon  tha  valociLy  ia  that  glva.n  by  Lha  aualed  l-KT-.STDCKFV) 
modal.  Tha  ahoak  valociLy  agrofu  wlLli  Lhu  LaaL  data  In  th.ia  ranga, 
which  aarvaiJ  auauntial l.y  an  a  conoiatancy  cliijck  baaauuii  tha  prodlcLad 
overprauaurua  wara  loniid  Lo  ’'grau  thara  with  tha  taut  data.  lha  uhock 
arrival  t.tma  lor  tha  acalad  1-KT-.STT)(HFV)  mode.''.  Ir  about  0.6  nil lliaucondu 
lasH  than  iiiuaunrud ,  which  la  axpactad  aa  it  la  a  nuclaar  bl.aat  wava 
mot' 0.1 . 

liutwean  tiia  2A.!i-f':  station  and  the  array  cunr.er  the  predicted 
shock  velocity  along  tlie  radial  increaues  becansa  ul  th.i  greater 
utrangth  0.I  Llio  M2  shock  that  forms  at  274.5  Loot.  Ilia  .(.LI'IjFCI-I  Mq 
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Hhock  1h  p,ouoL'ai.Ly  HLron^'r  Lima  mimaurod ,  buL  by  Lhu  Liiim  llu-  iaclclonL 
Hhock  in  auar  Llio  array  caaLar  Lho  aaparntJon  ia  arrival  Lima  baLwaan 
tlia  p  rad  i  a  Lad  tihock  and  Llia  LaaL  data  luui  dacrauaad  by  about  U.b  milli- 
Haaoadn . 

ilia  (ihoak  1  ruat  LraJaaLoriuH  aioa}',  tlia  blnaatorM  ara  aomparad 
la  Flppira  \2.  At  tlia  Hll-lt  iiLatioa  tlia  KKKLKCT-i  Hhaak  arrlvaa  aba/id  oi 
tlia  ti'HL  ahuaku  by  about  (>.2  mil  1  I aaaoadn .  iiia  aaparatioa  daarauaan 
Hiip.hLiy,  to  aliouL  O.U  mlillm'aoada  t't  about  tha  bO-it  atatJoii,  and  IL 
raiiiaiaa  auHaiitlaliy  aoiuUaat  irom  Lliara  on  Inward  to  ilia  array  caatar. 

It  la  aoLawortliy  Lliat  tlia  pradlatail  apaad  oi  tlu*  laaldaat  nlioak  iaarauai'ii 
la  tha  iiaal  law  iaat  iruiii  Llia  aaatar  roup, lily  aimilmr  Lo  tha  apaadup 
aliowa  by  tlia  taut  data, 


SECTION  ^ 
DISCUSSION 


Tho  ovi'i'in'CHHuro  mouHurjnxjiiLH  I'fom  the  MIGHTY  MACH  11-2 
(MB  11-2)  !-ilx-l)urnL  loHl:  luive  been  employed  here  nu  tlie  banin  I'or 
evaluation  ol  the  K1'',1''1,1''.C1T-')  code  reauitM.  llel'ore  revlewii.j;  the.  com- 
parlHon,  a  lew  eoinmenta  are  noted  regard  1  ni’  the  tiae  of  the  teat  re.sull!!. 

The  detonation  ul  nix  larpe  hi}>li“uxploHive  ehurtten  to  produce 
blant  waven  nuirtclently  identical  In  npact"  and  time  to  converjie  at  the 
array  center  with  the  at'curacy  that  mip.ht  occur  ior  nuclear  explosionn 
in  a  dillMcuit  technicai  tank.  Evidence  ol  the.  problem  can  he  ueen  in 
the  dlll'erencen  obtained  In  tlu'  overprennure  hintorien  muanured  on 
dll'ierent  blnectorn  tor  the  Interaction  ol  nimply  two  blant  waven, 
l''i),'uru  Va.  The  shock  overprennure  ior  the  pair  intercepting  on  the. 
.V-tO-dec,  binector  In  'Jli  percent  lower  than  Ior  the  pair  intercejit Inj>  on 
the  (lO-dep,  hlnectiH',  and  the  latter  in  25  and  30  iiercent  below  the  two 
vaiuen  coiiii)uted  iroiii  the  biant  models  (Hectlon  3-2), 

The  overprennure  tent  liicordn  are  prenented  In  l''i(',ui'en  (i  and  V 
p,roui)i'd  la  paliii  by  nlmllai'  Location!)  ior  c.omixu'iiion  pur'ponen.  The. 
dliierencen  between  tracen  ior  any  one  pair  in  readily  ai)pnrent.  Tin' 
diiierencen  aie  In  arrival  time,  waveiorm  and  maynituile, 

Many  oi  the  diiierencen  tliat  could  be  expected  between  the 
bliint  waven  themnelven  —  nucli  an  time  of  firrival,  nphericity,  etc. 
would  ('.eneraUy  be  expected  to  tend  to  lower  the  interaction  over- 
prennuren.  A  hydi'oc.ode  c.al.culation,  on  the  other  liund,  havinji  identical 
blast  waven  and  thiit  emi)loyed  accurate  hydrodynamic  modellinp,  --  uucli  an 
"nhai'p"  Kankine-llupiOn lot  niiock  waven,  snificiently  small  cell  or  mesh 
nl/,en,  etc.  —  would  be  expected  to  p.ive  overprensures  that  would  occur 
under  ideal  burnt  juid  atmonpheric  conditions,  A  jp)od  hydi'oc.ode  cal¬ 
culation  In  expec.t)!d  to  nerve  in  deiininp,  the  upper  limit  ot  over- 
pri'timires  tlnit  could  be.  expected  ior  blant  WJive  interact  J  onn ,  in 
partlculai'  for  nucie.'tr  burntn. 


The  Rl'U'’LHGT-'i  code  in  its  present  iorm  is  expected  to  be  most 
accurate  prior  to  interaction  ot  the  refracted  R[  waven  at  the  charge 
railiain.  Thin  period  terminates  with  the  interception  of  the  triple 
point  li-HfM^  and  the  charge  radial,  figure  Ac,  which  in  i)redicted  to 
occur  on  radliiin  at  2A.5  feet  from  the  array  center.  During  this  period 
the  shock  overpresnuren  predicted  along  a  charge  bisector  for  the 
Incident  shock  Mj^  are  1  ound  to  agree  well  with  the  maximum  overpressures 
mu.'inured  from  82  feet  inward  to  A1  feet  in  the  MB  II-2  tests.  Figure  8. 
The  2A.5-ft  intercept  of  the  triple  point  and  the  charge  radial  is  in 
good  agreement  with  the  2A  feet  predicted  by  Kuhl  (Reference  1)  and  the 
test  renults.  Definite  conclusions  regarding  the  overpressure  wave 
forms  predicted  for  this  period  cannot  be  made  because  of  some  anomalies 
in  the  tent  data  for  points  on  tl'.e  charge  bisectors,  Figures  7a-e. 
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However  there  Lh  good  agreenient  lietweee  tlie  predleted  wave  I’orm  and  one 
test  record,  for  Ulie  point  of  first  intercept  at  284  feet  on  the  60-deg 
eluirge  bisector,  iMgure  7a,  starting  at  one  miHisecond  after  shock 
arrival,  and  there  is  supportive  agreement  with  other  recau'ds. 

Tlie  KHKhKCT-d  prediction  of  tlte  R|-shocK  arrival  time  for 
points  on  the  charge  radial  is  In  good  agreement  witli  the  test  results 
for  the  82  and  41-Ct  stations,  Figures  6a-b.  Tlie  ItFFLECT-3  prediction 
has  the  shock  .iniiip  at  the  l)eg inning  of  the  R2  wave,  becanse  it  employs 
tile  sharp  R|  shock.  Tlie  test  records  show  a  very  sharp  pressure  spike 
for  this  R2  wave  following  the  shock  and  the  iiredicted  results  essentially 
fair  through  these  wave  forms.  Resolution  ol  the  overpressure  jicnk 
would  require  some  reduction  in  the  cell  siaes  used  in  the  code. 

The  R2  and  Ibj  shocks  (Flipire  4e)  have  not  lieen  modelled  in  tlie 
present  form  of  the  RKFLliCT-'l  code  as  sharp  Rankine-llugoniot  shocks,  so 
these  shocks  become  smeared  over  several  cells  in  the  results.  This 
smearing  causes  the  pressures  ahead  of  these,  shocks  to  rise  early.  For 
example,  the  rise  of  the  incident  M|  shock  overpressure  predicted  along 
a  charge,  bisector,  Figure  8,  commencing  at  the  22~ft  station  and  amounting 
to  about  10  percent  is  attributed  to  this  R2-slu)ck  snuuiring  effect.  The 
predicted  incident  shock  overpressure  is  found  to  be  higher  than  measured 
by  up  to  30  percent  irom  18  i/2  ieet  inward  to  lU  feet  which  is  attri¬ 
buted  in  part  to  the  smeared  R2  shock.  The  marked  ci. hub  in  the  M)  shock 
overpressure  coimiieneing  at  the  y-ft  stiiLlon  is  associated  with  the 
arrival  of  the  bulk  of  the  smeared  R2  shock. 

The  R2  shocks  intereeptlig  charge  bisector  form  the  R} 

waves.  The  early  arrival  of  the  leading  edge  of  the  smeared  R2  waves 
appears  to  affect  the  predicted  overpressure  wave  forms  at  stations  on 
the  bisector  from  82  feet  inward  to  10  feet,  Figure.  7c-g.  From  82  feet 
inward  to  42  feet  it  appears  to  cause  a  pressure  rise  about  0.01  to  0.02 
seconds  early,  beginning  at  about  0.00  seconds.  At  20  and  10  feet  the 
overpredietions  beginning  immediately  with  the  arrival  of  the  incident 
Ml  shock  are  attributed  to  it. 

After  the  predicted  formation  of  the  incident  M2  shock  (at 
24.5  feet)  on  the  radial  and  the  incident  M3  shock  (at  y  feet)  on  the 
bisector,  the  incident  shock  overpressure  predicted  by  the  REFLECT-3 
code  climbs  moaotonically  along  botli  the  radial  and  the  bisector  inward 
to  the  array  center.  Figure  8.  This  result  is  in  contrast  to  the  trend 
g.Lven  by  the  shock  interaction  analysis  of  Reference  1  of  decreasing 
shock  overpressures. 

Relative  to  this  question  the  test  results  show  the  following. 
The  predicted  Incident  shock  verpressure  at  22.5  feet  on  the  radial 
agrees  well  with  the  prediction  of  Reference  1  and  the  test  results  at 
20  feet.  Figure  8.  Comparisons  of  predicted  wave  forma  and  test  records 
have  not  been  made  for  the  radial  Inside  of  20  feet,  so  an  evalu.ation 
cannot  be  given.  At  10  feet  on  the  bisector  the  advance  arrival  of  the 
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smeared  shock  affects  the  comparison.  But  by  4  feet  this  effect 
would  be  expected  to  be  essentially  over,  and  at  4  feet  and  0  feet  the 
predicted  incident  shock  overpressure  agrees  well  with  the  test  results, 
Figure  7h-i.  It  is  concluded  from  this  comparison  that  tiie  strength  of 
tlie  incident  shock  increases  after  the  M2  and  M^j  shocks  are  formed,  but 
the  conclusion  is  tentative,  dependent  upon  calculations  with  sharp  R2 
and  tJbocks  and  smaller  culls. 


Tlie  r.iain  ref  leered  (M-R)  wave  formed  from  the  implosion  of  all 
six  blast  waves  at  the  array  center  is  of  primary  concern  as  it  results 
in  tile  highest  overpressures  observed,  Figures  9  and  10.  At  0.43  feet 
from  the  center  the  peak  overpressure  predicted  for  the  M-R  wave  by  the 
REFLECT-3  code  is  about  23  percent  lower  than  measured.  Figure  10.  At 
4  feet  on  the  bisec  .or  the  predicted  value  falls  slightly  below  the  data 
and  at  10  and  20  feet  on  the  bisector  and  radial,  respectively,  it  lies 
somewhaL  above.  'I'liere  is  thei'elore  a  tri'iul  for  l.he  prod  ( c.  ted  peak 
overpressure  of  tlie  M-R  wave  to  decrease  more  slowly  than  for  the  test 
results  and  to  he  higher  at  the  further  distances  than  measured  in  the 
test . 


The  test  records  indic.ute  a  rather  sharp  peak  for  the  M-R 
wave,  particularly  for  the  array  center.  The  predicted  M-R  wave  forms 
are  stretclied  out  more  in  time.  This  stretcli-out  is  expected  to  affect 
tlio  peak  value  predicted.  Tlie  stretch-out  is  attributed  to  the  finite 
ceil  sixe,  and  it  is  expected  to  have  a  greater  effect  on  the  comparison 
close  to  the  array  center  where  the  cell  dimension  becomes  more  important 
and  the  test  record  indicates  tlie  most  rapid  decay  timewise. 

A  reduction  in  the  cell  size  used  in  the  REFLECT-3  code  would 
tend  to  raise  the  peak  pressures  predicted  in  the  M-R  wave,  the  in¬ 
crease  could  well  result  in  raising  the  peak  pressures  of  the  M-R  wave 
at  all  distances  from  the  center.  In  fact,  the  predictions  might  well 
then  all  fall  well  above  the  measured  peak  overpressures.  Without 
further  information  on  the  characteristics  of  the  imploding  shock 
segments  in  a  multiple-burst  event  relative  to  the  effects  of  non- 
uniformity  in  blast  detonations,  atmospheric  and  terrain  effects,  etc., 
it  is  reasonable  to  speculate  that  pi!rhap.s  higlier  peak  overpressures 
might  be  attainable.  This  posslhility  Is  supported  by  the  rapid  decay 
indicated  in  the  measured  overpressure  in  Figures  9  and  10,  a  factor  of 
two  in  a  distance  of  only  4  to  6  feet.  Figures  9  and  10,  which  might 
imply  stringent  time-of-arrivui  reciuirv.ments . 

In  general  it  is  concluded  that  the  potential  of  the  REFhECT-3 
code  has  been  well  demonstrated  in  these  comparisons  with  the  MISERS 
BLUFF  II-2  data.  The  two  improvements  to  the  code  that  have  bee.n 
indicated  throughout  are  (1)  to  represent  the  reflected  shocks  R2  and 
Rg  as  sharp  Rankine-Hugoniot  shocks,  as  are  the  R]^  reflected  shock  and 
the  three  incident  Mach  shocks,  M-j ,  M2  and  M-^,  in  the  present  REFLECT-3 
code,  and  (2)  to  reduce  the  cell  sizes  for  better  wave  form  definition. 
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SECTION  5 
CONCLUSIONS 


'rlie  REFLEC'i'-3  code  has  been  developed  by  making  limited 
modifications  to  the  original  REFLECT  code  (DNA  Report  No.  3470F) .  The 
objective  has  been  to  determine  from  comparisons  with  the  test  results 
of  the  MISERS  BLUFF  II-2  event  the  potential  of  the  code,  if  modified, 
for  computing  the  interaction  of  multiple  nuclear  blast  waves  from 
ground-level  detonations.  The  MISERS  BLUFF  I1--2  test  consisted  of  the 
simultaneous  detonation  of  six  120-ton  ANFO  charges  placed  in  a  hexa¬ 
gonal  array  with  100-ft  separations.  The  following  conclusions  were 
reached  from  the  comparisons. 

1.  Between  100  and  41  feet  of  the  array  center  the  predicted 
incident  shock  overpressures  for  the  charge  bisector  agree 
well  with  the  maximum  values  measured. 

2.  Within  about  22  feet  of  the  array  center  tlie  smearing  of  the 
R2  and  R3  refracted  shocks  affects  the  predictions  of  the 
incident  shock  overpressure.  Overpredictions  relative  to  the 
test  results  of  up  to  30  percent  are  observed.  Within  4  feet 
of  the  array  center,  where  the  smearing  effect  may  be  less, 

the  predictions  agree  well  with  the  incident  shock  overpressures 
measured . 

3.  The  peak  overpressure  predicted  for  the  main  reflected  wave  is 
1615  psi  at  0.43  feet  from  the  array  center  (the  closest  point 
from  the  center  computed) ,  which  is  only  about  23  percent 
below  the  test  result.  At  4,  10  and  20  feet  from  the  center 
the  agreement  is  very  good.  Based  on  considerations  concerning 
the  code  calculations  and  test  factors,  it  is  possible  tiiat 
higher  overpressures  might  be  obtainable  than  measured;  it  is 
expected  that  they  could  be  adequately  predicted  following 
code  improvements  indicated  below. 

4.  The  two  improvements  to  the  code  that  have  been  indicated 
throughout  are  (i)  to  represent  the  two  reflected  shocks  R2 
and  R^  as  sharp  Rankine-Hugoniot  shocks  and  (ii)  to  reduce  the 
sizes  of  the  cells. 
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APPENDIX 

THE  PFAB  CODE 


GENERAL  DESCRIPTION  OF  CODE. 

The  PFAB  computer  code  was  developed  as  a  modification  of 
the  FAB-2D  code  of  Reference  9  for  computing  the  main  reflected  wave 
in  the  six-burst  problem.  The  method  employed  is  to  solve  the  planar 
two-dimensional  transient  fluid-flow  problem  of  the  reflection  of  an 
initially  prescribed  blast  flow  from  the  corner  of  a  wedge-shaped 
boundary,  as  indicated  in  Figure  13. 

The  flow  field  outside  of  the  computational  cell  system  shown 
in  Figure  13  is  a  prescribed  time  dependent  flow,  which  is  in  the  present 
application  given  by  the  REFLECT-3  code.  While  the  REFLECT-3  code  does 
compute  this  same  problem  as  part  of  its  late  time  solution,  the  PFAB 
code  is  more  realistic  in  that  the  main  reflected  blast  wave  is  repre¬ 
sented  as  a  sharp  Rankine-Hugoniot  shock  front  (which  is  smeared  in  the 
REFLECT-3  code)  and  the  PFAB  geometric  cell  layout  is  better  adapted 
to  the  features  of  the  main  reflected  wave. 

The  code  computes  the  flow  in  a  moving  cell  coordinate  system 
similar  to  the  one  shown  in  Figure  13.  The  cell  system  consists  of  an 
arbitrary  number  of  "sectors",  three  are  indicated  in  Figure  13  (designated 
J=1  to  3)  each  of  which  is  subdivided  into  an  arbitrary  number  of  cells, 
as  illustrat(jd  in  the  figure  for  sector  2  (J“2)  for  a  3  radial  cell  con¬ 
figuration  (1=1  to  5) .  Cell  dimensions  and  distances  are  expressed  in 
terms  of  the  radial  distance  r  and  the  polar  angle  0,  with  the  origin  of 
coordinates  being  taken  at  the  corner  of  the  wedge-shaped  space. 

The  outer  boundary  of  the  cell  system  represents  the  outer 
limit  of  the  calculated  flow  region,  whd.ch  is  initially  specified  and 
moves  radially  thereafter  for  each  sector  according  to  the  speed  of 
the  shock  wave  at  the  outer  boundary  of  the  outermost  cell  of  the  sector 
and  the  local  outer  prescribed  flow  region.  The  inner  and  outer 
boundaries  of  all  cells,  such  as  aa’ ,  bb'  in  Figure  13,  are  taken  to 
move  outward  similarly,  at  speeds  proportional  to  the  ratio  of  their 
initial  radial  distances  from  the  origin  of  coordinates  to  the  initial 
shock  front  radius  for  the  same  polar  angle. 

Within  each  sector  the  radius  from  the  origin  of  coordinates 
to  the  curved  circumferential  cell  boundaries  (aa',  etc)  is  taken  to 
vary  linearly  with  the  polar  angle  between  the  values  at  the  two  rays 
designating  the  sector  boundaries. 


These  shapes  are  not  precisely  true  sectors  in  that  their  outer 
boundaries  deviate  slightly  from  circular  arcs. 
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The  pressure,  density  and  velocity  components  in  each  cell 
(abb'a'a,  etc.) »  as  expressed  in  polar  coordinates,  are  taken  to  be 
constant  throughout  the  cell. 

The  flow  between  adjacent  cells  in  the  moving  coordinate  system 
is  computed  using  the  Godunov  technique  in  the  same  manner  as  for  the 
REFLECT-1  code  (Reference  2).  More  specifically,  in  the  interior  of  the 
disturbed  region,  the  flux  conditions  at  cell  boundaries  are  computed 
as  locally  isentropic  shock-expansion  phenomena  for  large  pressure  and 
velocity  differences  between  cells,  using  an  isentropic  exponent 
(Y  =  (dln(p) /dln(p) )  )  which  is  the  average  of  the  values  for  the  two 

adjacent  cells;  for  small  pressure  and  velocity  differences  between  adja¬ 
cent  cells  a  linearized  approximation  is  used.  For  the  outermost  cell 
boundary  of  the  disturbed  region,  the  shock  front  velocity  and  associated 
fluxes  are  calculated  on  the  basis  of  the  exact  Hugoniot  relationships 
for  air. 

The  code  assumes  a  real  air  medium  with  thermodynamic  equation- 
of-state  properties  given  by  Brode's  1965  analytical  representation 
(Reference  4) . 

Input  to  the  code  consists  of  the  specification  of  the 
pressure,  density  and  velocity  in  the  entire  disturbed  fluid  field  at 
a  time  shortly  before  the  incident  blast  wave  strikes  the  corner  and 
also  the  time  dependent  incident  blast  flow  conditions  outside  of  the 
computed  flow  region  for  all  later  times.  These  required  initial  and 
time  dependent  flow  conditions  would  be  obtained  from  the  results  of  a 
REFLECT-3  run.  The  conditions  would  be  inputted  to  the  PFAB  code  by 
an  interpolation  code  designated  REFRA  (Reference  10) .  The  REFRA  mating 
code  was  developed  for  use  with  the  REFLECT-1  code.  Minor  changes  are 
required  for  its  use  with  the  REFLECT-3  code. 

Program  output  consists  of  tables  of  pressure,  density  and 
velocity  and  configuration  geometry  printouts  for  selected  times. 

The  PFAB  code  has  been  set  up  on  a  CDC  CYBER  176  computer. 

The  code  has  a  core  requirement  of  164K  octal  SCM  and  376K  octal  LCM 
for  a  5000  cell  configuration. 
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